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We describe a method for recovering of the substructure, internal dy- 
namics and geometrical shapes of clusters of galaxies. Applying the method 
VO . to the Virgo cluster, we first, reveal the substructure of the central 4 arc 

^^ I degree field of the Virgo cluster by means of S-tree technique. The existence 

f^ ■ of three main subgroups of galaxies is revealed and their dynamical char- 

acteristics are estimated. Then, using the previously suggested technique 
(Ref. Q), the bulk flow velocities of the subgroups are evaluated based on 
1-^ ■ the distribution of the redshifts of the galaxies. The results enable us also 

to obtain a secure indication of the elongation of the Virgo cluster and its 
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O I positional inclination. 
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X ■ 1 Introduction 



The present paper aims to find out a method enabling the study of the 
subgroupings and their bulk flows, as well as the shapes of the clusters of 
galaxies, based on rather general assumptions. The existence of subgroups 
of clusters in general is a known fact, revealed by different methods for var- 
ious clusters (see e.g.a). The study of a sample of Abell clusters from ESO 
Key program (ENACS) performed by S-tree method enabled to conclude the 
existence of subgroups - galaxy associations, dynamical entities with remark- 
able properties as a possible common feature of clusters of galaxies 13 'O'H. 
We use the method to analyze the substructure and internal dynamics of 
the central region of Virgo cluster. Though the Virgo cluster because of its 



close location is one of the best studied clusters, the precise measurements 
of distance have been performed for relatively few galaxies which have been 
used for the accurate determination of the Hubble constant. These observa- 
tional studies at present are being intensely continued (see e.g. Ref . [^] , [0] ) , 
so that the most accurate study of the internal dynamics of the cluster is 
especially desirable, since each galaxy can participate a bulk flow motion, 
apart from the Hubble expansion. 

To study the internal dynamics one has first to analyze the hierarchical 
structure of the cluster. By this first step we reveal the substructure of 
the Virgo cluster core within the ±2 arc degree field centered on M87 using 
the S-tree method. On this point the present study is the generalization of 
the previous studyO, where the Virgo cluster within central 1 arc degree 
field has been studied. For the second step - the reconstruction of the bulk 
regular flow of the subgroups, we use the method developed in Ref.|j|] based 
on the data on the radial velocities of the members of the cluster. This 
method has been already applied to reveal the bulk motions within the Local 
Group of galaxies, and hence to obtain the apex and the 3D velocity vector 
of the motion of the Local Group with respect to the Cosmic Microwave 
Background frame H. 

For the Virgo cluster our analysis indicates the existence of 3 main sub- 
groups of galaxies. Then, we estimate the transversal components of regular 
motions of the subgroups with respect to each other. Among the main results 
of this study we have to mention the unambiguous indication of the elon- 
gated nature of the Virgo cluster, and combined Owith the rest kinematic 
parameters of the system we could even get constraints on the inclination 
angle and the degree of elongation. 

The ongoing studies on better estimation of the distances of individ- 
ual galaxies, particularly using the Cepheids, Tully-Fisher methods (e.g. 



iiJ'LLU'LL^ , the method of Surface Brightness Fluctuations for Virgo's bright- 
est galaxies Ii[, and other methods will enable from the results obtained here 
to determine more precisely the 3D bulk flow of each subgroups, and thus 
the mean motion of the Virgo cluster itself. Essential complementary in- 
formation on the parameters of the clusters and groups of galaxies can be 



provided by the X-ray datali^ . 



2 Virgo cluster substructures 

To study the substructures of the Virgo cluster core we used the data from 
the CfA redshift catalogueJiS. We dealt with the four arc degree central 
field within: N4316, RA= 12h 20m 12s and N4584, RA= 12h 35m 46.8s with 
the Virgo center coordinates RA=12h 27m 50.4s, DEC=12d 55m 55.2s. The 
sample includes 146 galaxies with available redshifts and V magnitudes. 

For t h e description of S-tree method one can refer not only to the original 
works'tL^'tliSl, but also to the recent applications of this method to various 
clustersliJl'O, where it is discussed in details. Therefore here we only stress 
that this method is revealing the degree of mutual interaction of members 
of the system using self-consistently the information both on the 2D coor- 
dinates and redshifts of individual galaxies of the cluster as well as on their 
magnitudes. 

The results of our study of the Virgo cluster core showed the existence 
of three main subgroups of galaxies containing 54, 17 and 37 galaxies, cor- 
respondingly. The lists of those galaxies along with their coordinates and 
redshifts are given in Tables 1, 2 and 3. 

While comparing the results of this analysis with those of the 1 degree 
fieldlEl it is interesting to notice, that we see the strong correspondence with 
the 1 degree field membership results. This fact is remarkable since it shows 
the robust character of the results of subgrouping derived by S-tree method 
even when limited areas of the clusters are studied. Relative stability of the 
subgrouping at various scales follows also from physical considerations. 

As seen from Table 4 below, though the mean redshifts of the subgroups 
are rather different, their redshift dispersions are quite close to each other. 
While the study of the regular motions of the subgroups will be done in 
next sections, already the remarkable difference of line-of-sight velocities is 
already indicating their possible significant radial bulk velocity with respect 
to each other. 

3 Bulk flow reconstruction 

The basic idea of the reconstruction procedure of the 3D motions within an 
N-body system from data on ID line-of-sight velocities relies on the existence 
of correlation between the velocities of different members of an interacting 
system. From this point of view our scheme of 3D velocity reconstruction 
and the S-tree method have identical physical background. For the case of 



Table 1: VIRGO subgroup I: 54 galaxies. 



Name 


I 


b 


^CMB 


Name 


I 


b 


^CMB 


N4316 


280.722 


70.963 


1599 


N4371 


279.666 


73.369 


1273 


N4374 


278.183 


74.475 


1363 


N4377 


275.310 


76.206 


1700 


N4379 


273.742 


76.967 


1394 


N4380 


281.900 


71.820 


1298 


N4390 


281.819 


72.270 


1452 


13328 


282.325 


71.897 


1337 


13331 


280.458 


73.575 


1608 


A1223+1308 


279.110 


74.545 


1594 


N4411 


283.890 


70.818 


1616 


13344 


278.454 


75.247 


1702 


N4411 


284.091 


70.856 


1605 


N4417 


283.452 


71.520 


1178 


13356 


281.390 


73.389 


1433 


13371 


282.534 


72.783 


1260 


N4429 


282.361 


73.012 


1463 


N4431 


280.995 


74.132 


1243 


13374 


283.584 


71.970 


1202 


N4436 


281.176 


74.175 


1454 


A1225+900 


285.316 


70.800 


1442 


N4451 


285.083 


71.334 


1195 


A1226+94 


285.139 


71.512 


1374 


N4459 


280.123 


75.842 


1540 


13413 


283.586 


73.469 


1735 


N4477 


281.538 


75.610 


1680 


N4478 


283.377 


74.392 


1698 


N4479 


281.855 


75.576 


1183 


N4483 


286.771 


71.223 


1215 


A1228+1219 


284.137 


74.159 


1578 


N4486 


283.761 


74.489 


1620 


N4497 


285.185 


73.806 


1427 


13457 


284.374 


74.818 


1796 


N4503 


286.081 


73.413 


1688 


13468 


287.006 


72.516 


1703 


13468 


285.632 


74.122 


1438 


A1230+926 


288.084 


71.497 


1563 


N4515 


280.529 


78.319 


1258 


N4516 


283.225 


76.739 


1280 


13487 


288.443 


71.748 


1489 


N4519 


289.177 


71.049 


1562 


13499 


287.655 


73.337 


1555 


N4528 


287.628 


73.675 


1702 


13510 


287.986 


73.446 


1703 


13518 


289.268 


72.045 


1771 


13520 


285.897 


75.816 


1416 


N4540 


283.618 


77.801 


1605 


A1232+927 


289.947 


71.644 


1617 


N4551 


288.150 


74.681 


1523 


N4564 


289.541 


73.920 


1492 


13583 


288.272 


75.709 


1448 


13602 


291.863 


72.660 


1607 



Table 2: VIRGO subgroup II: 17 galaxies. 



Name 


I 


h 


^CMB 


Name 


/ 


h 


^CMB 


N4321 


271.133 


76.899 


1884 


N4330 


278.752 


72.906 


1885 


N4383 


272.099 


77.758 


2015 


N4405 


273.413 


77.585 


2063 


13349 


280.227 


74.227 


1801 


N4421 


275.764 


77.021 


1925 


13369 


274.990 


77.549 


2048 


13392 


278.247 


76.760 


2001 


1796 


276.345 


78.126 


1919 


N4474 


280.797 


76.005 


1949 


N4486B 


283.392 


74.563 


1914 


13457 


284.374 


74.818 


1796 


13470 


286.240 


73.508 


1829 


13501 


285.437 


75.594 


1932 


13586 


289.098 


75.005 


1892 


N4578 


291.682 


72.125 


2613 


N4579 


290.374 


74.361 


1845 











Table 3: VIRGO subgroup III: 37 galaxies. 



Name 


/ 


h 


^CMB 


Name 


/ 


h 


^CMB 


13239 


278.239 


73.229 


979 


N4328 


271.511 


76.943 


841 


N4387 


278.828 


74.465 


913 


A1223+1513 


275.513 


76.439 


830 


N4402 


278.760 


74.783 


565 


13363 


280.318 


74.349 


1120 


N4424 


283.865 


71.389 


767 


A1224+936 


284.169 


71.332 


1047 


N4435 


280.148 


74.889 


1101 


N4440 


281.363 


74.170 


1068 


N4442 


284.140 


71.819 


849 


N4445 


284.625 


71.482 


635 


13381 


282.244 


73.722 


967 


13393 


281.313 


74.817 


794 


N4452 


282.696 


73.729 


553 


N4458 


281.078 


75.151 


1011 


13412 


284.978 


72.086 


1098 


A1226+1243 


282.451 


74.438 


866 


N4469 


286.087 


70.905 


833 


N4486A 


283.984 


74.386 


778 


N4491 


284.808 


73.636 


826 


13459 


284.941 


74.359 


606 


13466 


285.451 


74.028 


1114 


N4506 


283.826 


75.578 


1006 


1798 


281.397 


77.480 


734 


N4523 


283.106 


77.362 


582 


13517 


289.588 


71.588 


764 


13522 


284.036 


77.476 


980 


N4548 


285.669 


76.830 


871 


13540 


287.350 


75.326 


1057 


N4550 


288.087 


74.635 


706 


A1233+1239 


288.036 


74.793 


608 


N4552 


287.914 


74.967 


647 


A1233+1408 


286.827 


76.240 


1078 


13578 


289.978 


73.592 


1021 


N4571 


287.486 


76.659 


663 


A1235+1508 


288.207 


77.363 


1055 











stellar systems (large N) the problem of reconstruction of 3D velocity distrf 



bution of stars has been solved decades ago by Ambartsumiantll , without 
any a priori assumption on the form of the distribution function. The only 
natural assumption was its translation invariance, i.e. the phase space dis- 
tribution of the system can be split as 

dP = ^s{xi, X2, X3)dxidX2dx3 X ^K{vi,V2,V3)dvidV2dV3, (1) 

where ^sl^^i) 2:^25 2^3) is the 3D spatial distribution of objects and <1>k(wi5 ^2, ^'s) 
is the 3D velocity distribution function. The efficiency of such reconstruc- 
tion procedure is however closely related to the number of objects considered. 
For the case of galaxy clusters (small A^), it turns out that Ambartsumian's 
method cannot be applied directly. However some quantities of interest, 
for example the bulk flow of individual clusters, can be extracted from the 
data by assuming reasonable form of the velocity distribution function $k • 
Herein we assume that the velocities distribution of galaxies inside the clus- 
ter can be described by gaussian random isotropic components of velocity 
dispersion a^ plus the 3D mean peculiar velocity of the cluster, i.e. 

^k(^^i, ''^2, -i^a) = 9{vi; Vi,av) g{v2] K, a^) g{v3; V3, a^), (2) 

where Vr is the radial component of the cluster bulk flow and Vi and V3 
are its components perpendicular to the line-of-sight. Note however that 
this hypothesis does not imply any severe limitations in our problem. More- 
over it is indeed justified by the fact that due to exponential instability of 
gravitating systems and their mixing properties, the correlations in physical 
parameters of particles have to split. This effect can be followed in numerical 
experiments. 

Another limitation comes from the fact that galaxies inside a cluster 
take part in the Hubble flow, implying then that a fraction of the observed 
redshift is due to the spatial elongation of the cluster along the line-of- 
sight. A common practice is to assume that this effect is negligible, i.e. all 
the galaxies lie at the same distance. In Ref.[Q] we used a more physical 
assumption allowing a 3D spatial extension of the cluster with gaussian 
isotropic distribution around the cluster center, i.e. 

$s(a;i,a;2,X3) = g{xi;0,as) g{x2;0,as) g{x3;0,as), (3) 

where as is the spatial dispersion of the cluster. Under this working as- 
sumption it is shown in Appendix A3 that the observed probability density 



reads in terms of the redshift z and the angular position with respect to the 
center of the cluster = {6i, 63) of the galaxies as 

dPohs = giz-{z);Vi Oi + Vs 63, dobs) X 5(^1; 0, ag,) giOs; 0, ae^,) dO^dOidz (4) 

where a'l^^ = o"^ + <7|, {z) is the mean redshift of the cluster and (Tq^ and gq^ 
are the angular sizes of the cluster in the plane perpendicular to the mean 
line-of-sight. It thus turns out from Eq. @ that the tangential velocity of 
the cluster VJ;an = (^1,^3) can be evaluated by using standard statistical 
technique. 

4 Application to Virgo subgroups 

For each Virgo subgroup the system of coordinates {xi,X2-,X3) is defined as 
follows. The X2-axis is directed towards the mean angular position of sub- 
group galaxies {{l),{b)), the xi-axis is parallel to the galactic plane and the 
X3-axis is chosen in a way that the {xi,X2,X3) system forms a direct trihe- 
dron. The xi, X2 and X3 coordinates then are deduced by transforming the 
cartesian galactic coordinates {z cos / cos b,z sin / cos b,z sin b) so that the two 
trihedrons coincide. Because the velocity dispersion a^ is small compared 
to the mean redshift (z) of the cluster, one has within the approximation of 
small angles xi « 61 (z) and X3 ^ 63 (z) where 61 and 63 are the position 
angles along the xi and xa-axis respectively, and X2 ~ z. 

The 3D data distribution of Virgo subgroup I is given in Figure 1 , where 
the X2 coordinate has been transformed to X2 = 2; — {z) for convenience. 
Redshifts are expressed in km s^^ in the CMB frame. The influence of 
peculiar velocities is clearly seen on X1X2 and X22;3 projections. Indeed 
the data distribution is much more elongated along the X2-axis due to the 
presence of the velocity dispersion. The two left panels of the figure reveal 
a correlation between 9i and 63 (or equivalently between xi and X3). Such 
a correlation which cannot be described by Eq. (0) clearly means that we 
have to modify our initial hypothesis, namely that the spatial distribution 
of the subgroup galaxies admits a central symmetry. 

FIGURE 1. 

Indeed the hypothesis of a central symmetry is stronger than it is nec- 
essary for the problem we are dealing with i.e. for the tangential bulk flow 
estimation. In Appendix A3 it is shown that the sufficient condition for 



Table 4: Estimate of the parameters for the 3 Virgo subgroups. Units: (l) 
and (6) in degrees, crg^, ag^ and a'g in radians and (z), Uobs, ^21 and ^423 in 
km s~^ 





Virgo I 


Virgo II 


Virgo III 


A^gal 


54 


17 


37 


(Z) 


1497.91 


1958.22 


862.96 


(l) 


284.162 


281.446 


283.635 


(b) 


73.658 


75.637 


74.522 


(Te^ 


0.0182 


0.0302 


0.0176 


^ez 


0.0326 


0.0338 


0.0331 


-e. 


0.0165 


0.0183 


0.0173 


Pl3 


-0.172 


-0.631 


-0.033 


Ai3 


-0.23 ±0.07 


-0.71 ±0.14 


-0.10 ±0.09 


fobs 


167.22 


171.26 


178.95 


A2I 


2485 ± 1375 


2684 ± 2264 


-489 ± 1703 


A23 


938 ± 767 


476 ± 2023 


-109 ± 904 



evaluating Vtan is the absence of correlation between the spatial distribution 
of galaxies along the line-of-sight and perpendicular to it, i.e. 021 and 023 
are vanishing in Eq. (|46|). In this case, the observed distribution function 
takes the form 



dPobs = 9{z- {z);A2i9i+A23 03,(^ohs) 

xg{ei;Ai3 e-s, o-^J 5(6*3; 0, ag.^) de-idOidz 



(5) 



where the tangential velocities are Vi = A21 and V2 = A23. Here the pa- 
rameter Ai3 allows us to include into consideration the correlation between 
6*1 and ^3, as required by the observational data. 

The estimates of tangential bulk flows of the three Virgo subgroups are 
given in Table 4 using the formulae derived in Appendix A3. Due to the 
relatively small number of galaxies in the subgroups the error boxes are 
large, but nevertheless are informative. Note, that for subgroups I and II 
the bulk flow peculiar tangential velocities are anomaly high. These values 
are significant i.e. are not due to the sampling errors. 



5 The spatial structure of Virgo subgroups 

The statistically significant correlation between 6i and 9^ for Virgo sub- 
groups I and II advocates against the hypothesis of central symmetry for 
the spatial distribution of their galaxies. This means that these systems rep- 
resent spatially elongated configurations. Their elongation can be described 
by the following ellipsoidal 3D spatial distribution function 

^s(yi, 2/2,2/3) = 9 {yr, 0,ai) g{y 2; 0,a2) g{y 3; 0, as) (6) 

where at least one of the dispersions, say ai, differs from the two others. 
In this case, the orientation of the proper system of coordinates (2/1,^/2,2/3) 
of the ellipsoidal distribution has to be defined with respect to the frame 
of analysis, namely (xi,X2,X3). This is done by introducing 3 rotation an- 
gles Q, /3 and 7, where the angles a and (3 are respectively the longitude 
and the latitude of the yi-axis in the (xi,X2,X3) frame and 7 is the angle 
between the intersection of the planes X1X2 and y22/3 'ind the y2-axis. As it 
is shown in Appendix A2 the spatial probability density in the coordinate 
frame (xi,X2,X3) has to form 

dPs = g{x3;0, ai) g{xi; ai3 X3, an) x g{x2;a2ixi + a23X3,a{i^) dxiX2X3 (7) 

where the parameters 013, 021, 023, o"!, an and a'm are linked to the proper 
spatial dispersions o"i, (T2 and (T3 and the orientation angles a, (3 and 7 by 
the formulae given in that appendix. The case 021 = and 023 = arises if 
the ellipsoid is symmetrical with respect to the plane perpendicular to the 
line-of-sight (a = and 7 = 0). 

A it is shown in the appendix A3, the observed distribution function has 
the same form as in Eq. (g|), but then the parameters A21 and ^23 are 

A21 = Vi + 021 Hq D ■ A23 = V3 + a23 Ho D (8) 

where Hq is the Hubble's constant and D is the distance of the considered 
subgroup. These expressions show us that in general it is not possible to 
disentangle between the presence of tangential subgroup bulk flow and the 
effect of inclination along the line-of-sight of the ellipsoidal distribution of 
subgroup galaxies. In particular, it is very likely that the large values of A21 
and A23 obtained for the Virgo I and II subgroups are partly due to this 
inclination effect. 

Finally, let us consider the spatial structure of Virgo su bgroups I and 
II based on the results given in Table 4. The initial model contains 10 free 



parameters, namely the 3 bulk flow components Vr, Vi and V3, the velocity 
dispersion a^, the 3 proper spatial dispersions cJi, (T2 and (T3 and the 3 
orientation angles a, [5 and 7. Our approach enables us to estimate the 6 
independent parameters ^13, ^23, ^21, o"e , (^83 and dobs- Assuming that the 
subgroup is at rest in the CMB frame, i.e. 14 = 0, Vi = and V3 = 0, one 
can deduce for a given value of the velocity dispersion ay the characteristics 
of the spatial structure of the subgroup. 

We have performed this experiment for Virgo subgroups I and II based 
on the data in Table 4. Figures 2 and 3 show respectively the variation of 



parameters ui, (T2, CT3, a, fi and 7 as a function of ctjjj = J cr'^^^^ — o"^ for the 
subgroups I and II. Since the problem does not admit a general analytical 
solution, these curves have been obtained numerically. 

The results show that both subgroups reveal an elongated structure with 
significant difference between ui, 02 and us, whatever is the contribution of 
the true velocity dispersion a^ in the observed redshift dispersion Uobs- 

FIGURES 2 and 3. 

6 Conclusion 

Thus the technique we had developed for 3D velocity reconstruction of 
galaxy configurations H notwithstanding to the inevitable error boxes, can 
be successfully applied to the Virgo cluster subgroups. 

By a first step we have estimated the tangential bulk flow velocity of the 3 
subgroups under the assumption that the observed redshift dispersion is not 
due to their 3D spatial structure. For two of the subgroups, namely Virgo I 
and Virgo II, these values were found statistically significant but anormaly 
high (2656 km s~^ and 2726 km s~^ respectively). This fact indicates the 
validity of our main working hypothesis, i.e. that the clustered galaxies 
exhibit an elongated spatial distribution along the line-of-sight. 

We have shown that the knowledge of galaxies redshifts alone still does 
not permit to disentangle between the presence of tangential bulk flow and 
the effect of inclination of an elongated 3D spatial structure. Nevertheless, 
while both effects are certainly contributing to our derived statistics, we can 
assert unambiguously that both subgroups are spatially elongated, with the 
specific characteristics of their 3D structures depending on the kinematic 
properties of the subgroups. 
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A Useful formulae 

The probability density of a random variable x given as a gaussian distri- 
bution with a mean xq and dispersion ax 

dPc = g{x;xo,crx)dx = —^= — exp — ^ — ] dx (9) 

V2TTax V 2^^ / 

fulfills the properties 

g{x + A; xq, a^) dx = g{x; xq - A, a^) dx (10) 

g{Xx; Xxo,\X\ ax) d{\X\ x) = g{x;xo,ax) dx, (11) 

where A is a real number. The following relation also holds 

g{x; xi,ai) g{x; X2,a2) = g{x; xq, ctq) g{xi;x2, a), (12) 

with xo, Co and a defined as 

2 2 

O- =0-1 +0-2 , ^0 = 2 , 2 • y^^) 

a I i- ^2 



Xo = 2 , 2 ^1 + 2 I -2 ^2- (14) 



o-f + a| erf + fj^ 

We describe an ellipsoidal 2-dimensional distribution in the system of 
coordinates (yi,j/2) by the probability density 



dPeii = 5(2/1; 0,cri)dyi X g{y2;0,a2)dy2 (15) 
or using the properties (p!0|), ( [llD and (^) in the system of coordinates 

(2;i,X2) 

dPeii = g{xi; Ax2,a'i) g{x2;0,a2) dxidx2. (16) 
with a'2, a'l and A defined as follows 

/2 2 2,2-2 / '^l'^2 /i„N 

(T2 = 0-2 COS a -|- o"]^ sm a ; o"i = — ^ (17) 
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A = sin a cos a | ^ — • (1^) 



\ ^2 / 

Here the angle a characterizes the rotation Ra transforming (xi,X2) into 
(yi)2/2) coordinates, i.e. 

_ J yi = cosa xi + sina X2 .,„s 

I 2/2 = — sm a xi + cos a X2. 

B The spatial distribution 

The 3-dimensional spatial distribution of considered galaxies can be ex- 
pressed in the system of coordinates (2/1,2/2,1/3) as 

dPs = 5(2/1; 0, cTi) dyi X 51(2/2; 0, 0-2) dy2 x ^-(2/3; 0, 0-3) dys. (20) 

The spatial orientation of this trihedron is herein characterized by 3 rotation 
angles a, (3 and 7 with respect to the frame of analysis (xi,X2,X3). The 
angles a and (3 are respectively the longitude and the latitude of the 2/1-axis 
such that coordinate transformations read as 

2/1 = cos a xi + sin a X2 
P-a '■ { y2 = —sina xi + cosa X2 (21) 

2/3 = ^3, 

where (2/1 > 2/2; 2/3) i^ the new system of coordinates after a rotation R^ of 
angle a around the X3-axis, 



Rf- 



y[ = cos P 2/1 + sin P 2/3 

2/2 = 2/2 (22) 

2/3 = - sin /3 2/" + cos /3 2/3 



Here (2/1,2/2)2/3) is the system of coordinates after rotating on an angle (3 
with repect the 2/2 axis, and the angle 7 is defined by a rotation R^ around 
y'l axis so that 2/2 and 2/2 axes coincide, i.e. 

2/1 = 2/1 
R-^i : { 2/2 = cos 7 2/2 + sin 7 2/3 (23) 

2/3 = -sin7 2/2 + cos 7 2/3. 
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Our aim is to express the density probability of Eq. ( poD in terms of the 

coordinates of analysis, namely {xi,X2,X3). Using (|l6|) and in view of (|25 

one has 

dFi = g{y2;0,a2)dy2 g (y 3; 0,(73) dys 



= 5(2/2; -4 2/3, 0-2) 5(2/3; 0' ^3) dy'2dy'3, 
with (T3, (72 and A defined as 



/ 2 9 9 9 9 

'^3 = "^s cos 7 + (T2 sin 7 



/ 0"2<73 



A = sin 7 cos 7 



^2-^3 



cj: 



/2 



The relation ([l|) and definitions (|2|j2l|,|3|) yield 

dF2 = g{yi;0,ai)dyig{y'^;0,a'3)dy's 

= 5(2/1 ; B X3, a'{) g{x3] 0, 0-3) dy'{dx3, 



with (T3, o"^' and i? defined as 



af = a'^ cos^ /3 + a? sin^ /3 



<Ti 



(TifT3 



(24) 

(25) 
(26) 

(27) 



(28) 



B = sin /? cos (3 



r^r 



cj: 



//2 



(29) 



Using the relations ( [10| ) and (pLlD, the definitions (21,2^ and defining 61, 62, 
^3) ^4) ^5 and tg as 



61 = sin a ^2 = cos a + A sin a sin f3 

65 = — cos a 63 = sin a — A cos a sin P 
b^ = A cos/3 ^6 = B. 

we obtain for the following probability density 



dF3 = 5(2/2; ^ 2/3' 0-2) 5(2/1'; ^ X3, a'{) dy'2dy 



= gih X2; bs xi + 64 X3, CJ2) (i|62X2| 
X5((6i X2; 65 2;i + be 2;3, o"'/) c^l^sa^^il 

= 5(^1^2 a;2; bib3 xi + 6164 X3, |6i|o-^) d|6i62X2| 
xg{b2bi X2; 6265 a;i + 6266 X3, |62|o-i ) d\b2b5Xi\. 



(30) 
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Using dl^ ) and since 6163 — 62^5 = 1 the spatial distribution of Eq. (^) can 
be rewritten in terms of the coordinates {xi,X2,X3) 

dPs = g{x3;0,ai)g{xi;ai3X3,aii) 

(31) 

Xg{x2] 021 Xi + 023 X3, Cr(jj) dxiX2X3, 

where the parameters 013, 021, 023, o'l, an and cr^jj are obtained from Eqs. 



(pq25,25,2e,29) 



4 


= bW,' 


^hla'C 


5 


ai = 


II 

--03 


021 = 


0-ii2 


o"iii = 


0"II 


; «23 


1 
CTll2 



ai3 = ^2 fte - &i ^4 (32) 

(63^)2^^ - 65 &i4') (33) 

(6462^1^ - ^6 61 CT^') (34) 

The problem of extracting the values of parameters (a,/3,7,(Ti,(T2,o"3) from 
a given set of observed parameters (013,021, 023, a"i,<Tii,(Tiii) is not straightfor- 
ward in general. The special case of a 2-dimensional spatial distribution 
deserves mentioning however. If one of the dispersions in Eq. (|2^) vanishes 
(e.g. a\ = 0), the previous calculations lead to 



J' n . ,/2 ^11 . ,/2 01 
C^i — U , cr2 ~ . 2 ' '^3 ~ 2/3 

sm^a cos2/3 




COS [} 

A = (013 + cos Q tan /?) ; B = — tan/3, 

sin a 




which implies 




tana = ; tan/? = — 023 sina 

021 


(35) 


and finally yields 




1 / 22 n\ 

, ., / NN 1 / ,9 cm cos /) \ 
sinhlntan7 =-- A^ - 1 + ^ —^ 
A y ai^ sm^ a J 


(36) 


2 f. ^ \ ^i^ 

/T — 1 1 


(37) 


V tan 7/ cos^/) 


2 
ai = (l-yltan7)-^. 


(38) 
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C Distribution in the phase space 

The phase space distribution of the considered dynamical system can be 
spht as fohows 

dPps = dPK X dPs, (39) 

where dPs is the 3D spatial distribution of Eq. (pO) and the probability 
density dPK describes the distribution of peculiar velocities v = {vi,V2,V3). 
We chose the system of coordinates such that the 2;2-axis points toward the 
mean angular position of the cluster and such that the xi-axis is parallel to 
the galactic plane; vi, V2 and ^3 are thus components of the galaxies peculiar 
velocity respectively along xi, X2 and X3 axes. In this frame the 3D velocity 
distribution is given by the 3D mean peculiar velocity of the cluster plus 
random isotropic components of velocity dispersion a^, i.e. 

dPK = g{vi;Vi, ay)dvi g{v2;Vr, ay)dv2 g{v3;Vi, a^)dv3, (40) 

where Vr is the radial component of the bulk flow and Vi and V3 are its 
components perpendicular to the line-of-sight. The distance of the cluster 
is D and the angular positions with respect to the center of the cluster are 
9i and ^3. In the approximation of small angles, 61 and 9^ are 



Then the redshift z of a galaxy expressed in km s^^ units is 
z = HqD + v.u 

7^ HqD + V2+ HqX2 + 0iVi+ 63 V3, 



(41) 



(42) 



where Hq is the Bubble's constant and u is the line-of-sight velocity. 

Integrating the probability density (iPx over the two unobserved tangen- 
tial components of the velocity vi and V3 we have 

dP, = g{z; HoD + Vr + HoX2 + Oi vi + 63 V3,ay) dz (43) 



Where ([10|), ( [TlD and (|T2) have been used. It follows from this result and 
from the expression Eq. (3^) for dPs the integral dPps in Eq. (|39| ) over the 
unobserved variable X2 yields 

dPohs = g{z; HoD + Vr + A21 9i + A23 O3, CJobs) 

(44) 
xg{ei;Ai3 63, a'gjg{e3;0, ae^) de3deidz, 
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where the parameters A13, ^421, A23, o"obs, (^e^ and cjg^ are related to 013, 
^21 5 0231 CI"!, ail and a'm from Eqs. (|32| , ^^ ) . This characterize the spatial 
structm'e of the cluster, and velocity dispersion a^ and cluster's tangential 
velocities Vi and V3 through the following formulae 

'^obs = ^v^ + ^'m ; ^13 = ai3 (45) 



""ei = TTT; ' ^^^3 = TT-F^- (47) 



A21 = Vi + aai Ho D ; 7123 = ^3 + 023 ^0 D (46) 

So far as the angles 9i and ^3 have been defined in a way that their averages 
over the sample vanish (i.e. (^1) = and (^3) = 0), the observed probability 
density dPohs of Eq. (E3) can be rewritten as 



O^^obs = 9{Z; A21 01 + A23 6*3, CTobs), 

(48) 
xg{ei-Ar3 63, a'g^ ) 5(^3; 0, ae, ) dO^ddidZ, 

where the observable Z is defined as 

Z = z-{z) = z-{HqD + Vr). (49) 



From Eq. ( |4q ) it turns out that the parameters ^13, A21 and A23 can be 
obtained using a standard multiple regression technique, i.e. 

Coy {01,63) ,_„, 

Cov(6l3,6'3) 

_ Cov(gi, g3)Cov(gi, Z) - Cov(gi, gi)Cov(g3, Z) 

'' Cov(ei,03)Cov(^l,e3)-Cov(^l,ei)Cov(03,^3), ^ ' 

_ Cov{63, 93)CoY{ei,Z) - Cov{6i,93)CoY{e3, Z) 

'' COV(^1,01)COV(^3,^3)-COV(01,^3)COV(01,^3), ^ ^ 

where Cov(x, y) = {{x — {x)){y — {y))) is the covariance of random variables 
X and y. Estimates of the dispersions dobs; f^e-i and a'g^ are obtained using 

ae,^ = 00^(63,93), (53) 

4i' = ^ei' -2Ai3 Cov{9i,03) + A?3 ^^32, (54) 

cTobs^ = fTz^ - 2^21 Cov(^i,Z) - 2 ^23 Gov (03, Z), (55) 
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+Al^ cj0,2 + Al^ ae^^ -2 A21 ^23 Cov(03, ^i) 
, where gq^ and az are defined as 

ae,2 = Cov(^i,0i) ; a/ = Cov(Z,Z). (56) 

Finally, the standard deviations i.e. accuracies of the estimators of parame- 



ters Ai3, ^421 and ^423 given Eqs. ( p0|j5l|j5^) can be represented as a function 



of the number N of galaxies in the cluster 

A^i3 = 4^ -^, (57) 

yN 0-6)3 

AA23 = -^ , ^ ^^, (58) 



AA21 = ^= ^^ — , (59) 



ll3 = 


1 ^ei 
^ ^63' 




1 


1 


O"obs 


Viv 


Vi - Pli 


f^es 


1 


1 


O-Qbs 



where /331 is the correlation coefficient between variables 9i and ^3, i.e. 

pL = Cov(^.^3)Cov(.„g3) c,,(,^^,^) C0V(^3, ^3) = ^?3 ^ (60) 



^^1 
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Figure captions. 

Figure 1. 

Virgo subgroup I: 54 galaxies. The top left panel visualizes the 3D in- 
formation contained in the data. The ^1^3 plane is perpendicular to the 
mean line-of-sight ((/) = 284.162,(6) = 73.658). The redshift depth is pro- 
portional to the radius of the circles, backside galaxies are shadow, frontside 
plain. The 3 other plots are the projection of this 3D distribution on plane 
X1X3, X2XS and 2;iX2 respectively where X2 = z — (z) . 

Figure 2. 

Virgo subgroup I: Variation in function of a[^^ of the parameters char- 
acterizing the ellipsoidal spatial structure (o"i, (72, cr^) and its orientation 
(a,/?, 7). The dispersion a'm reads in terms of Uobs = 167.22 km s^^ and the 

i^ _ ^2 

'v 



velocity dispersion a^ as a'm = A/o"^bs — cr, 

Figure 3. 

Same caption as fig. 2 but for Virgo subgroup II: The observed dispersion 
is o"obs = 171.26 km s~^. 
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